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ABSTRACT 



An important number of multiwavelength studies of low ionization nuclear emission-line regions (LINERs) were ded- 
icated to investigate the excitation mechanism responsible for the detected emission lines. Radiative emission from 
accretion into a super massive black hole (SMBH) is now believed to be, in an increasing number of LINERs, the 
source of excitation. However, the accretion mode is not yet firmly understood, and could be explained in either a 
thin accretion disk or a radiatively inefficient accretion flow (RIAF). Our purpose is to study the X-ray properties of 
LINER sources with definite detection of a broad Ha emission line in their optical spectra, LINER Is from Ho et al. 
sample. These objects preferentially harbor a low luminosity active nucleus at the center and show small or no intrinsic 
absorption (<10 22 cm _1 ). We compare their X-ray properties to both X-ray binaries and luminous AGN. We analyzed 
all available X-ray archived XMM-JVewton and Chandra observations of 13 LINER Is satisfying the above criterion in a 
systematic homogeneous way. We looked for any correlations between the X-ray properties and the intrinsic parameters 
of our sample of LINER Is. An absorbed power-law gave a good fit to the spectra of 9 out of the 13 sources in our 
sample. A combination of a thermal component and an absorbed power-law were required in the remaining 4 sources. 
We found a photon index for our sample between 1.3 ± 0.2 for the hardest source and 2.4lg'3 for the softest one with 
a mean value of 1.9 ± 0.2 and a dispersion a = 0.3. The thermal component had a mean temperature kT~ 0.6 keV. 
Significant short (hours to days) time-scale variability is not common in the present sample and was observed in only 2 
sources (NGC 3226 and NGC 4278). Three other sources indicate a possible variability with a low K-S test probability 
(2%-4%) that the nuclear emission originates from a constant source. On the other hand, significant variability on a 
longer time-scale (months to years) is detected in 7 out of the 9 sources observed more than once. No significant Fe Kq 
emission line at 6.4 keV was detected and upper limits were derived for the 4 sources with a high enough signal to 
noise ratio around 6 keV. Finally, we established, for the first time for a sample of LINER Is, that the photon index L 
is significantly anticorrelated to L2-10 kev/LEdd- Whereas this anticorrelation is in contrast to the positive correlation 
found for type 1 AGN, it is similar to the one seen in XRBs in their low/hard state where a radiatively inefficient 
accretion flow is thought to be responsible for the X-ray emitted energy. 

Key words. Accretion, accretion disks - galaxies: active - galaxies: nuclei - X-rays: galaxies 



1. Introduction 



It is widely known that a large fraction of local galaxies 
contain emission line nuclei that are the result of low nu- 



length studies were attributed t o this subject, loo king for a 



radio, sometimes variable, core (Nagar et al.|2005) or a vari 



able UV core (Maoz et al. 2005 



clear activity. Ho et al.] ( 1995 1997a|b ) showed that more 
than 40% of 486 nearby galaxies with Bt < 12.5 mag 



oy gaiaxies wirn r>T 
could be considered as active with optical spectra classi- 
fied as Seyfert nuclei, low-ionization nuclear emission-line 



regions (LINERs, Heckman 1980 ), or transition objects (ob- 
jects having intermediate spectra between LINERs and H II 
nuclei). Pure LINER sources would be the most common 
component representing 20% of all 486 galaxies. The ion- 
ization mechanism responsible for the excitation of emis- 
sion lines in LINER sources is an ongoing matter of de- 



in nearby LINER sources. 
Nevertheless, the most used tool to search for an active nu- 
cleus in a LINER is to look for a hard 2-10 keV unresolved 
core that could not be due to diffuse em ission from shock 
heated gas or from unusually hot stars (|Terashima et al 
[20001 iHoet al.||2001|[Dudik et al.l|2005] |Flohic et al.||20Qi , 
|Gonzdlez-Martfn et al.|2006||2009[7Zhang et al.|2009| ). How 
do LINERs harboring a low luminosity active nucleus com- 
pare to luminous Seyfert galaxies and quasars? 



(Dopita & Sutherland 


1995|, starburst activity (Alonso- 


Herrero et al.|2000 


Terievich & Melnick|1985 


) , or a low iu- 



minosity active Galactic nucleus (AGN)l Many multiwave 



Maoz] (2007 1, using high angular resolution multiwave- 
length observations of 13 LINER sources, demonstrated 
that the luminosity ratios in different wavebands, mainly 
UV to X-ray and radio to UV luminosities, follow the same 
trend as luminous Seyfert galaxies. The authors did not 
find any sharp change in the spectral energy distribution 
(SED) of their sample of 13 LINERs compared to more 
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luminous Seyfert and quasar nuclei, suggesting that a thin 
accretion disk ma y persist at low accretion rates. Moreover, 
Pian et al. ( 2010[ ) detected up to 30% flux variations on half 
a day time-scale in 2 (NGC 3998 and M81) out of 4 LINER 
and low luminosity AGN sources observed in X-ray with the 
XRT onboard Swift. They combined their X-ray fluxes with 
simultaneous UV fluxes coming from the UVOT instrument 
and showed that the SED and the UV to X-ray flux ratios 
of their 4 sources sample are consistent with those of more 
luminous sources and that LINERs may have similar ac- 
cretion and radiative processes at their center compared to 
luminous Seyfert nuclei. 

On the other hand, the faintness of LINER sources com- 
pared to luminous Seyfert galaxies and quasars has been at- 
tributed to a different accretion mechanism owing to some 
observational contrast between the two classes. No broad 
nor narrow Fe Ka emission line at 6.4 keV have been de- 
tected in the spectra of the LINER sources with the highes t 
signal to noise ratio ( |Ptak et al.||2004j |Binder et aL"l|2009| ), 
X-ray short ti me-scale v ariability has been detected in only 
a few sources ( Ptak et al.|1998 Awaki et al.|2 001), and the 
disappearance of the bi g blue bump in the UV band in the 
SED of LINER sources |Ho|1999[[2008| , all these signatures 
could indicate the disappearance of the thin accretion disk 
at low luminosities and that a different accretion mech- 
anism is responsible for the emission in LINER sources. 
It has been suggested that when the mass accretion rate 
falls below a critical value M C rit, the density of the disk 
could become too low for radiative cooling to be effective. 
The trapped heat will expand the thin accretion disk into a 
pressure-su pported thick disk with a very low radiative e ffi- 
ciency (see |Quataert|200lj|Narayan fe McClintock|2008[ for 
reviews). Such radiativeiy inefficient accretion flow (RIAF) 
models successfully explained the spectra l energy distribu- 
tion of a large n umber of LINER sources ( Ptak et al.|2004| 
Nemmen etaT1|2006| |2010[ ). 

Another way to assess the geometry of the accretion 
mode in AGN is to compare them to their less massive 
counterparts, X-ray binaries (XRBs) . Shemmer et al. ( 2006 ) 
showed that the X-ray spectral slope, T, of Seyfert 1 galax- 
ies and quasars and the Eddington ratio, iboi/^Bddi ar e 
positively cor related, similar to X RBs in their high/soft 

and references therein) 



2008, 



state (see also Shemmer et al. 
Such a behavior could be explained in an AGN accretion 
mode consiste nt with an optically thick geometrically thin 
accretion disk (|Shakura fe Sunyaev|1973[).|Gu fe Cao|(|2009| 



i9) 
lid 



performed a similar study on a broad sample of LINERs an< 
low luminosity Seyfert galaxies. They found a significant 
anticorrelation between F and the Eddington ratio for the 
local Seyfert galaxies in their sample analogous to XRBs in 
the low/hard state where a RIAF mode of accretion takes 
place. However, no strong correlation was found when con- 
sidering only the LINER sources in their sample, owing, as 
suggested by the authors, to heterogeneous fitting models 
as they have collected their data from different studies. In a 



separate study, Constantin et al. ( 2009 ) analyzed the X-ray 
emission of a sample of 107 nearby galaxies including low 
luminosity Seyferts, LINERs, transitions (nuclei with spec- 
tra between Seyferts and LINERs), H II regions, and pas- 
sive galaxies (lacking optical emission-line activity), none 
of which show broad-line components. Using a Spearman- 
rank correlation, the authors found an anticorrelation for 
their sample between L and the Li, i/LEdd- By considering 
each class separately, a spearman-rank test showed that the 



anticorrelation persists for the different objects, except for 
the low luminosity Seyfert galaxies. 

Finally, broad optical emission lines, a characteristic 
property of classical Seyferts and quasars, are also found 
in nuclei of much lower luminosities. Thirty three sources 
out of the 221 n uclei classi f ied as Seyfert, LINER, or tran- 
sition objects in Ho et al. ( 1995 ) sample of nearby galax- 
ies show definite detection of a broad Ha emission, 16 of 
those (-17% of the total pure LINER sources) are LINERs 

LINER Is here- 



Ho et al. 1997b 



(noted as LINER 1.9 in 
inafter). In this paper, we are aiming for the study of the 
X-ray characteristics of these LINER Is observed with the 
current generation of X-ray telescopes, XMM-Newton and 
Chandra. Such a sample insures the responsibility of ac- 
cretion into a SMBH for the formation of the broad emis- 
sion lines (given the early type class of this sample where 
outflows form massive stars and/or supernovae are not ex- 
pected to be relevant), guarantees the non-existence of large 
obscuration, and enables X-ray comparison of this class 
with both XRBs and type 1 AGN. We introduce our sam- 
ple in § 2, § 3 represents the observations and the data 
reduction. Temporal and spectral results are given in § 4. 
In § 5 we discuss the results in the context of LINER ls- 
Seyfert-XRB connections, and a conclusion summarizing 
the main results is given in § 6. We report, in appendix A, 
some notes on the individual sources and in Appendix B 
we give spectral results to the surrounding sources around 
the centers of galaxies observed with Chandra. In the re- 
mainder of this paper, luminosities are calculated using the 
distances given in Table [I] derived with a Hubble constant 
H Q = 75 km s- 1 Mpc" 1 . 



2. The sample 

We select ed objects classifi ed as LINER 1.9 (LINER Is) 



sources mlHoetaL (1997b) showing a definite detection of 
a broad Ha emission line. This implies the definite existence 
of an AGN at the center of all of the sixteen selected galax- 
ies and its responsibility for the excitation of the detected 
optical emission lines. 

We excluded three sources from the sample: NGC 3642, 
NGC 4636 and NGC 1052. NGC 3642 did not have any 
archived XMM-Newton or Chandra observations. As for 
NGC 4636, all of the X-ray archived observations were stud- 



ied in extreme detail (Jones et al. 2002 



Ohto et al.|[2003l |(J'Sullivan et al-POOl^ 



2; 

B 



Xu et al, 
aldi ct al 



2002 



2009 



Xu et al.||2010 ) and show a complicated spectrum that re 



quires, for a good spectral parameters measurement, de- 
tailed imaging analysis and should be modeled including: 
sophisticated shock models, temperature and density gra- 
dients, and last but not least, steep abundance gradients in 
the core. Finally, the broad Ha emission line detected in the 
spectrum of NGC 1052 is attributed to polarization due to 
electro n scattering within the opening cone of an obscuring 
torus ( Barth et al.| 1999 ). NGC 1052 was classified as an 
obscured AGIN showing a large int rinsic absorption in the 



X-ray spectrum (Ng 
T999|. 



10 



23 



Guainazzi & Antonelli 



Tabic [T] shows the list of galaxies along with their corre- 
sponding right ascension and declination, Hubble type, dis- 



tance (taken from Tonry et al. [2001 otherwise from Tully 
1988), the mass of the black hole derived from the M-cr re- 



lation ( |Tremaine et al.|2002 Graham et al.|2010 ) where the 
velocity dispersion is taken from Ho et al. (20091. Multiple 
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Table 1. Properties of the 13 LINER Is showing definite detection of a broad Ha emission (taken from Ho et al. |l997b 
sample) . 



Galaxy Name 




R.A. 




Dec. 




nubble lype 


Distance 
(Mpc) 


log(M Bff ) 6 


NGC-266 


00 


49 


47 


8 


+32 


16 


40 


SB(rs)ab 


62.4 


8 


44 


NGC-315 


00 


57 


48 


9 


+30 


21 


09 


E+: 


65.8 


9 


06 


NGC-2681 


08 


53 


32 


7 


+51 


18 


49 


(R')SAB(rs)0/a 


17.2 


6 


78 


NGC-2787 


09 


19 


18 


5 


+69 


12 


12 


SB(r)0+ 


7.48 


8 


15 


NGC-3226 


10 


23 


27 





+19 


53 


55 


E2:pec 


23.6 


8 


06 


NGC-3718 


11 


32 


34 


8 


+53 


04 


05 


SB(s)a pec 


17.0 


7 


61 


NGC 3998 


11 


57 


56 


1 


+55 


27 


13 


SA(r)0? 


14.1 


9 


07 


NGC 4143 


12 


09 


36 





+42 


32 


03 


SAB(s)0 


15.9 


8 


18 


NGC-4203 


12 


15 


05 





+33 


11 


50 


SAB0-: 


15.1 


7 


73 


NGC-4278 


12 


20 


06 


8 


+29 


16 


51 


E1+ 


16.1 


8 


72 


NGC-4750 


12 


50 


07 


2 


+72 


52 


28 


(R)SA(rs)ab 


26.1 


7 


27 


NGC-4772 


12 


53 


29 


1 


+02 


10 


06 


SA(s)a 


16.3 


7 


46 


NGC-5005 


13 


10 


56 


2 


+37 


03 


33 


SAB(rs)bc 


21.3 


7 


79 



Note s. "Distan ces adapted from Tonry et al. ( 2001 1, otherwise from Tully (19881. Black hole mass cal culated using 



et al. (20101 updated M-a relation of Tremame et al. (20021 with stellar velocity dispersion taken from Ho et al. ( 2009) 



Graham 



snapshot observations (exposure time < 5 ks) were ex- 
cluded from the analysis either because of high background 
contamination (NGC 4143, obs.ID: 0150010201), low num- 
ber of counts detected (NGC 2787, obs. ID: 388), or severe 
pile-up (NGC 4203 and NGC 4278, obs. IDs:397 and 398, 
respectively). The final sample consists of 13 LINER Is with 
a total of 31 observations summarized in Table [2j 



3. X-ray observations and data reduction 

3.1. Chandra observations 

All of the LINER Is in our sample have at least one 
Chandra observation. Snapshot observations with an ex- 
posure time less than 5 ks were performed for eight 
sources (NGC 266, NGC 315, NGC 3226, NGC 3718, 
NGC 4143, NGC 4750, NGC 4772, and NGC 5005). Seven 
sources have observations with a sufficient exposure time 
for a detailed temporal and spectral study (NGC 315, 
NGC 2681, NGC 2787, NGC 3226, NGC 3998, NGC 4203, 
and NGC 4278). All of the Chandra observation s were ob- 
tained wit h the spectroscopic array (ACIS-S; Weisskopf 
et al.|2002 ) where the nucleus was placed on the aim point, 
except for NGC 3226, of the ACIS-S3 back-illuminated 
chip. They were taken in either Faint or Very Faint mode 
to increase their sensitivity. All of the observations are 
Chandra archival data obtained from chaser 1 . The log of 
the Chandra observations are listed in Table 2. 

All Chandra observations were reduced and analyzed 



version 3.175 2 (Broos et al.|2010). We started by using the 

~b 



in a systematic, homogeneous way (as in Younes et al 
2010 hereinafter Y10) using the CIAO software package 
version 4.2, Chandra Calibration Database, CALDB, ver- 
sion 4.3.1, and the ACIS Extract (AE) software package 



level 1 event file produced by the Chandra X-ray Center 
(CXC) to suppress the position randomization applied by 
the CXC Standard Data Processing when creating a level 
2 event file. We also corrected for the effects of charge- 
transfer inefficiency on event energies and grades. We fil- 
tered for bad event grades (only ASCA grades 0, 2, 3, 4 
and 6 are accepted) and hot columns to take account of 
several potential issues such as cosmic rays and bad pixels. 
Good time intervals, supplied by the pipeline, were applied 
to the final products. 

The LINER nucleus source position is determined af- 
ter running a wavelet transform detection algorithm, the 
wavdetect program w ithin the CIAO data analysis system 
( Freeman et al.|2002 ). This position is then given to the AE 
software that refines it, extract source photons, construct 
local backgrounds, extract source, and background spectra, 
compute redistribution matrix files (RMFs) and auxiliary 
response files ( ARFs) , by spawning the mkarf and mkacis- 
rmf routines of CIAO, and perform spectral grouping and 
fitting. 

Source events are extracted around the source centroid, 
inside a polygonal shape of the local PSF, generated by 
MARX 3 at the energy 1.497 keV using the ae_maie_psf tool 
implemented in AE. Background region is defined following 
the AE procedure. The background region is an annular re- 
gion centered on the source position where the inner radius 
delimit 1.1 x 99% encircled energy radius and the outer ra- 
dius is set such that the background includes between 100 
counts and 200 counts (depending on the brightness of the 
source). This background is obtained from a special image 



The 
age and 



ACIS 
User's 



Extract 
Guide 



software pack- 
available at 



http: / / cda.harvard.edu / chaser/Descriptions 



http://www.astro.psu.edu/xray/acis/acis_analysis.html 
6 http://space.mit.edu/ASC/MARX/ 
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Table 2. Log of the Chandra and XMM-Newton X-ray observations. 



Source Name 


Satellite 


Instrument 


Start Date 


Obs. ID 


Net Exposure-time 


l\ - Y i f ^ rice 

lNUU-2oo 


Chandra 


AUlb-b 


2001 June 01 


i c i n 

1610 


o n 


1NUU-315 


Chandra 


AUlb-b 


onnn /" ^ . . ,-, ,, no 

zUUU Uctober Uo 


000 


4. / 




Chandra 


AUlb-b 


zUUo February zz 


415o 


r r n 

55.0 




XJMJM-Jvewton 




onnc t , , l no 

2005 July 02 


nonconnom 

0305290201 


1 o o /oo o /oo o^ 

lo.z/zo.z/^o.z 


f~* OCOI 

lNUU-2o81 


Chandra 


AUlb-b 


2001 January 30 


one n 


on n 

80.9 




Chandra 


AUlb-b 


2001 May 02 


one i a 


7n n 
/9.0 


1NUU-2/8/ 


Chandra 


AUlb-b 


onn/i A.f~, r i o 

2004 May to 


A £ o n 


on o 
30.8 




AJMJM-Jvewton 


LLiu 


onn A I \ ■ ] -i f t 

2004 Uctober ID 


nonnocm ni 

UzUUzoUlUl 


i A o /oo cr /oo nfr 

14.2/28.5/28.0 


M/"! /"i OOO/' 1 

lNUU-322o 


Chandra 


AUlb-b 


i n nn I A i . ,. _ on 

1999 December 30 


8o0 


AG C 




Chandra 


AUlb-b 


2001 march 23 


1515 


o o 
2.2 




AJMJM-Jvewton 




onno A/r^«^x^ on 

2002 March 29 


ni ni n/inoni 
U1U1U4UoU1 


on c /o^ n /o^ n 
30.0/30.9/30.9 




AJMJM-Jvewton 


LLiu 


onno i ., nn 

2008 January 09 


n a nno^n i ni 

0400270101 


n A ofo 

94.3 


1NUU-3718 


Chandra 


AUlb-b 


onno T7''. . i . .. . . „ „ no 

2003 February Do 


onno 

3993 


t n 

4.9 




AJMJM-Jvewton 


LL1U 


onn/l a J„ no 

zl)U4 May (Jz 


nonn a oncrm 

0200430501 


n 

9.0 




AJMJM-Jvewton 


LLiu 


onn a at „ „ , i_. n /i 

zl)U4 iNovember U4 


nonn /t o i om <^ 

0200431301 


o o 

8.8 


NGL 3998 


Chandra 


AUlb-b 


onnc? ni 

zUUo July (Jl 


6781 


13.0 




AJMJM-lvewton 


LLiu 


onm A /r~ nn 

zUUl May 1)9 


nnnnnnni ni 

(juyuu^uiui 


O ft /if) CT /lO CT 

8.9/lz.5/12.5 


\in /^i A A AO 

1NUU 4143 


Chandra 


AUlb-b 


zUUl March zo 


151 ( 


o c 
Z.O 




A JM JM-J v e wton 




onno i\ t . , ... v \ , oo 

zUUo INovember zz 


ni £nm ncm a 
U15UU1U5U1 


no/ii n/ii n 

y.o/ii.y/ii.y 


1NUU-4203 


Chandra 


AUlb-b 


onnn in 

2UU9 March 1U 


i n c; o 


/ii n 

41.5 


IM"*" 1 /" 1 /I 0^70 

1NUU-42 i 8 


Chandra 


AUlb-b 


onnc T7. , l . ^. ....... . no 

zUUo February Uz 




O I .0 




Chandra 


AUlb-b 


onnr 1 A /r V. i 

2UUo March lb 


71)77 


lino 

11U.J 




Chandra 


AUlb-b 


onnc t,,k. oc 

zUUb July ZO 




C 1 A 

51.4 




Chandra 


AUlb-b 


onnc f \ . , v. , , ^. o yi 

zUUo Uctober 24 




i n r. i 
1U5.1 




Chandra 


AUlb-b 


onn'? TV 1 ,. i . . . ... on 

zL)U7 February zU 


f-7f\ O 1 

7081 


i i n t 

110.7 




Chandra 


ACIS-S 


2007 April 20 


7080 


55.8 




XMM-Newton 


EPIC 


2004 May 23 


205010101 


30.3/35.2/35.2 


NGC-4750 


Chandra 


ACIS-S 


2003 August 27 


4020 


4.9 


NGC-4772 


Chandra 


ACIS-S 


2003 February 14 


3999 a 


4.7 


NGC-5005 


Chandra 


ACIS-S 


2003 August 19 


4021 


4.9 




XMM-Newton 


EPIC 


2002 December 12 


0110930501 


8.7/13.1/13.1 



Notes. "Observations reported for the first time for the LINER 1 nucleus study. ^Exposure time corrected for solar flare intervals. 



where all events within the ~ 1.1 x 99% PSF circles of all 
the sources in the field were excluded (swiss cheese image). 
Background was modeled for snapshot observations. 

Piled-up observations were accounted for by excluding 
the core of the PSF (see Y10 for more details). We use 
the tool dmextract, called by the AE software, to create 
spectra over the energy range 0.5-8 keV. We used the tool 
ae-groupspectrum implemented in AE to group the spec- 
tra. Channels between 0.5 and 8 keV are grouped to have 
a three sigma (3er) signal to noise ratio, which corresponds 



to a minimum of 20 counts per bin, to enable the use of 
the x 2 statistics in the spectral analysis. The cash statistic 
(C-stat) is used to derive spectral parameters for snapshot 
Chandra observations with the background being modeled 
with t he cp linear background model developed by |Broos| 



et al. (2010). The background model is arbitrarily chosen 
to consist of continuous piecewise-linear (cplinear) func- 
tions with 2 to 10 vertexes. The model has 2 to 10 parame- 
ters representing the X-ray fluxes at the different vertexes. 
These vertexes are placed on the energy scale so that they 
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divide the energy range into intervals with approximately 
equal numbers of observed counts in the background spec- 
trum (0.1 to 10 keV). Vertex energies are chosen to coincide 
with the energies of actual events in the background which 
helps to prevent the vertex flux from being driven to the 
hard limit of zero during the fitting process (see § 7.5 of 
Broos et al.||20To| . 



3.2. XMM- Newton observations 

The log of the XMM-Newton observations is listed in 
Table 2. Eight sources were observed at least once 
with XMM-Newton (NGC 315, NGC 2787, NGC 3226, 
NGC 3718, NGC 3998, NGC 4143, NGC 4278, and 
NGC 5005) and two have multiple observations (NGC 3226 
and NGC 3718). In all of the o bservations, the EP IC-pn 
( Striider et aLpOOlJ and MOS flTurner et al||2001| ) cam- 
eras were operated in Imaging, Prime Full Frame or Large 
Window Mode (except for the long NGC 3226 observa- 
tion, where MOS cameras were operating in Small Window 
Mode 4 ) using the thin or medium filter. The Reflection 
Grating Spectra show only few counts for all the different 
observations and therefore they were not included in our 
analysis. We did not make use of the optical/UV dat a taken 
with the o ptical/UV monitor (OM) instrument (Mason 
et al. 2001 ) since this paper concentrates on the X-ray char 



acteristics of this sample. A multiwavelength study of our 
sample will be treated in a forthcoming paper. All data 
products were obtained from the XMM-Newton Science 
Archive (XSA) 5 and reduced using the Science Analysis 
System (SAS) version 9.0. Data are selected using event 
patterns 0-4 and 0-12 for pn and MOS, respectively, dur- 
ing only good X-ray events ("FLAG=0"). None of the EPIC 
observations were affected by pile-up, although severe inter- 
vals of enhanced solar activity, where the background count 
rate even exceeds the source count rate, were present dur- 
ing several observations (NGC 315, NGC2787, NGC 3226). 
In these cases, we reduced the background to 5% by ex- 
cluding the high background intervals which reduces the 
observation time usable for spectral analysis, sometimes to 
less than 30% of the raw exposure time (NGC 315). 

XMM-Newton source events of all of the LINER Is in 
our sample were extracted from a circle centered at the nu- 
cleus using two different radii of 10" and 25" 6 . We compared 
light curves and spectra of both extraction regions to check 
if any of the sources (jet emission, diffuse emission, and/or 
unresolved point-like sources) detected in the Chandra im- 
age between 10" and 25" around the nucleus contaminate 
the XMM-Newton nucleus emission (see Appendix B and 
online Fig. 1-4). No change is seen in the light curves, and 
the fit parameters of the two extracted spectra were con- 
sistent, within the error bars. Therefore, and to achieve 
better statistics and better constrain fit parameters, source 
events of all of the LINER Is in our sample observed with 
XMM-Newton were taken from the 25"-radius circle cen- 
tered on the nucleus. We added the spectral contribution of 
the different sources (jet emission, diffuse emission, and/or 

4 NGC 3226 is off axis during the long observation and hence 
is not detected with the MOS cameras when operating in a Small 
Window Mode 

http:/ /xmm.esac.esa.int/xsa/index.shtml 

6 The lower 10" limit was taken so to encircle at least 50% of 
the EPIC XMM-Newton PSF. 



unresolved point-like sources), derived from the Chandra 
observation and detected in a 25"-radius circle around the 
nucleus (see Appendix B for more details), to the spec- 
tral model used to fit the XMM-Newton spectrum. The 
particular case of NGC 4278 is discussed in detail in Y10. 
Background events are extracted from a source-free circle 
with a radius twice of the source on the same CCD. We 
generated response matrix files using the SAS task rmfgen, 
while ancillary response files were generated using the SAS 
task arfgen. The EPIC spectra were created in the energy 
range 0.5-10 keV to enable flux and model-parameter com- 
parison with Chandra spectra. They were grouped to have 
a signal to noise ratio of 3 with a minimum of 20 counts 
per bin to allow the use of the % 2 statistic. 



4. Results 

4.1. Light curves and hardness ratios 

Temporal analysis was only done for long exposure ob- 
servations, not including snapshot Chandra observations. 
Light curves and corresponding hardness ratios, defined as 
HR = (H - S)/(H + S), where S is the count rate in the 
soft 0.5-2 keV band and H is the count rate in the hard 
2-10 keV band, were extracted for all of the long obser- 
vations. We corrected the net count rate of the piled-up 
sources for the excluded fraction of the PSF. Chandra and 
XMM-Newton light curves were all binned with a time bin 
size of 1 ks for a reliable rms variability analysis. 

We first conducted a Kolmogorov-Smirnov, K-S, test to 
examine any potential variability within each observation. 
Based on this test, we do not find short time-scale (hours 
to days) variability in 5/10 sources (NGC 315, NGC 2681, 
NGC 3718, NGC 3998, and NGC 5005) with a K-S test 
probability >10% that the nuclear emission originates from 
a constant source. Three sources (NGC 2787, NGC 4143, 
and NGC 4203) indicate a possible variability with a K- 
S test probability between 4% and 2% that the nuclear 
emission originates from a constant source. Two XMM- 
Newton observations of two different sources exhibit sig- 
nificant short time-scale variability, both already reported 
in the literature, NGC 4278 (obs ID: 205010101, Y 10) and 
NGC 3226 (obs.ID: 0400270101, |Binder et al.]2009| , where 
the K-S test gives a probability less than 1% that the core 
emission is originating from a constant source. NGC 4278 
shows an increase at the beginning of the observation of 
10% on a time-scale of ~1.5 hour, the emission remains 
constant for the rest of the observation following that hint 
of variability. As for NGC 3226, variability is clear through 
the whole observation where a total increase of ~60% is 
detected bet ween the beginning and the end of the ~100 ks 
observation ( |Binder et al.|2009| . XMM-Newton light curve 



5] and the XMM-Newton and 
ler sources are given in online 



of NGC 3226 is shown in Fig 
Chandra light curves of the ot 
Fig. § and Fig. [7] 

To check more accurately any intrinsic variability am- 
plitude from the differe nt sources, we calcul ated the nor- 
malized excess variance ( |Nandra et aT|l997a[ ) for all of the 
long observations with the following expression: 



'NXS 



N^ 2 



N 

£ 

»=1 



[PQ- M ) 2 -a t 2 ] 



(1) 
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Fig. 5. Light curve (upper panel) and hardness ratio 
(lower panel) of the ~100 ks XMM-Newton observation of 
NGC 3226 binned to have a 1 ks resolution. The dashed 
lines show the averages on the count rate and hardness ra- 
tio. 

Table 3. Normalized excess variance for LINER Is with a 
relatively long exposure time. 



Galaxy Name 


2 

°NXS 


Exposure Time" 


°NXS 


Exposure time 




XMM-Newton 




Chandra 


NGC-315 


< 0.0027 


49.7 


< 0.0160 


55.0 


NGC-2681 






< 0.0400 


159.9 


NGC-2787 


< 0.0170 


37.1 


< 0.0453 


30.8 


NGC-3226 


(0.02±0.002) t 


94.3 


< 0.0700 


46.6 


NGC-3718 


< 0.0078 


18.4 






NGC 3998 


< 0.0001 


8.9 


< 0.0016 


13.6 


NGC 4143 


< 0.0060 


9.3 






NGC-4203 






< 0.0020 


41.6 


NGC-4278 


< 0.0012 


30.3 


< 0.0202 


470.8 


NGC-5005 


< 0.0042 


8.7 







Notes. "The exposure times in ks used to calculate the value 
of 

a Nxs- 6 Value corresponding to only the longest observa- 
tion of NGC 3226. Including the other XMM-Newton obser- 
vation would result in a cr% xs upper limit of 0.013. 



where N is the number of bins in a given light curve, 
Xi and <7i are the count rate and uncertainty of each bin, 
respectively, and \i is the arithmetic mean of the counting 
rates. To enable <J 2 NXS comparison between the different 
light curves, the bin size and the light curve segment du- 
ration should be taken equally. For that purpose, we first 
decided to use light curve segments of 20 ks, as usually 
done for luminous Seyfert galaxies, splitting any longer 
observations into multiple ones. That limited our sample 
to 6 sources observed with Chandra and 5 observed with 
XMM-Newton, not enough to draw any safe conclusions. 
Therefore, and owing to the heterogeneous sampling of the 
observations for this type of study, we decided to use the 
whole corrected exposure time of all of the long observa- 
tions. The mean of the <J NXS is taken for every source with 
multiple Chandra or XMM-Newton observations. The time 
bin size choice of 1 ks for all of the observations was taken 
to have a good signal to noise ratio with at least 20 counts 
in each bin and an acceptable number of bins in each light 
curve. 

Estimating the error on the uf^ xs could be a tricky 
task. The variability in an AGN light curve depends, on one 



10" 



t> 10" 



10 



NGC 3226t 



10 



10 



10 



10 



Fig. 8. (J xxs derived from the XMM-Newton observations 
as a function of the BH mass for our sample of LINER Is. 
Arrows represent upper limits. NGC 3226 is the only source 
showing clear short time-scale (~1 day) variability and thus 
a non- upper limit value on o J2 <[xs . 



hand, on the measurement errors of the data (e.g. Poisson 
noise) and, on the other hand, on the stochastic nature of 
the process underlyi ng AGN variability (e.g. red noise, see 
Vaughan et al.|2003 for a detailed discussion on this issue); 
even if a source is not intrinsically variable the mean and 
the variance of different light curves based on observations 
performed at different times will not be identical. We esti- 



mated error due to Poisson noise using the |Vaughan et aL 
(2003) equation 



err I a 



NXS) 




+ 



2a P 



N 



The unce rtainty owing to the r ed noise process has been 
presented by Vaughan et al. ( 2003 ) to depend on the power- 
spectrum s hape of the source wh ich we do not know a priori. 
Therefore, O'Neill et al. ( 2005[ ) estimated the error on the 
red noise process directly fro m the data. Our obs ervations 
are not well sampled to use O'Neill et al. (2005) method 
to determine any error due to the stochastic nature of the 
AGN X-ray variability. The error on <r NXS shown in equa- 
tion [2] was also used to estimate upper limits to the excess 
variance in the case of non-variability detection whenever 
a NXS ^ s negative or consistent with zero. Only one ob- 
ject, NGC 3226, in our sample shows a clear short time- 
scale variability during the longest 100 ks observation with 
a 2 NXS = 0.02 ± 0.002 (comp arable to the value of 0.014 
found by Binder et al. 20091). Upper limits were obtained 



for the rest of the sample. Table [3] an d Fig. [8] summarizes 
the results that are discussed in § |5.1[ 

4,2. X-ray spectral results 



The s pectral analysis was performed using XSPEC ( Arnaud 
1996) version 12.6.0. The photo-electric c r oss se ctions and 
the solar abundances of Wilms et al. (20001 are used 
throughout to account for absorption by neutral gas. An 
absorbed Galac tic column density de rived for every sin- 
gle source from Kalberla et al. (20051) (obtained with the 
W3NH tool 7 ) was applied to the different spectral models. 



http://heasarc.gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pl 
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Spectral uncertainties are given using A% 2 of 2.71, corre- 
sponding to 90% confidence for one interesting parameter, 
and to 95% confidence for upper limits. 



4.2.1. snapshot observations 

We began our spectral analysis with the study of the 
Chandra snapshot observations. Table [4] gives the best 
fit parameters to the snapshot observations of our sam- 
ple of LINER Is. Complicated models, like partial cover- 
ing and/or two power-law components could not be tested 
due to the low number of counts. A power-law modified by 
Galactic absorption gave a good fit in the case of NGC 266 
and NGC 4143. An additional intrinsic neutral absorp- 
tion was needed in the remainder of the cases. The pho- 
ton indices vary between 1.3 ± 0.2 for the hardest spec- 
tra to 2.1 ± 0.7 for the softest ones with a mean value 
of about 1.7. The hydrogen column density of the intrin- 
sic absorber had an upper limit of ~ 4 x 10 21 cm~ 2 in 
the case of NGC 4750 and NGC 5005. A value consistent 
with 10 21 < N H < 10 22 cm -2 was derived for the rest 
of the snapshot observations. In one case, NGC 3718, we 
find a somewhat larger column density with Njj in the or- 
der of ~ 10 22 cm' 



For NGC 50 05, a thermal component 
using the abundance ta- 



1985 



(meka l model, Mewe et al 

ble of |Wilms et al.||2000| with a 0.8+^ keV temperature, 
was included m the model to take account for some low en- 
ergy features, most likely due to diffuse emission from hot 
gas. In order to rigorously confirm the validity of our best 
fit spectral-parameter values derived using the C-stat and 
a modeled background, we compared our results to the re- 
sults derived from fits applying the x 2 statistics to all of the 
snapshots observations. Spectral parameters derived using 
the C-stat were all in agreement, within the error bars, with 
the results derived using the x 2 statistics; with smaller de- 
viation from the central value of one interesting parameter. 
We decided, as a consequence, to use the C-stat fits to cal- 
culate model fluxes in the soft, 0.5-2 keV, and in the hard, 
2-10 keV, bands. Table|6]gives the corresponding 0.5-2 keV 
and 2-10 keV observed fluxes and corrected luminosities. 






^ 0305290201: PN 






MOS1 






, MOS2 

k 








jJiililuiilijiJi 






mmf m° 

I I I 4156 



1 2 

Energy (keV) 






0400270101: PN ' ^ 1 \^ ^ ^^^^^^^ 






0101040301: p|n ' 




' ++.* tadllkkkitilfcj ... 


0101040301: MOS1 


'It t|'+tt1-H 


+ t \l l l l< |lll|l'l| l l | !l||l|i > l'l l l l l'l < l 




0101040301: JlO^ I 











Energy (keV) Energy (keV) 

Fig. 9. Upper left panel. Data and best fit model of the 
different spectra of NGC 315. A hardening in the Chandra 
ACIS spectrum (in blue) above ~1.5 keV is seen, relative to 
the XMM-Newton spectra (burgundy, green, and orange). 
Upper right panel. Residuals of the best fit model in terms 
of sigma. Lower left panel. Data and best fit model of the 
different spectra of NGC 3226. It is clear that more absorp- 
tion below 2 keV from cold material is taking place between 
the two XMM-JVewton observations (black representing the 
long ~100 ks observation). Lower right panel. Residuals of 
the best fit model in terms of sigma. 



4.2.2. Long-exposure observations 

We then turned to the analysis of observations with rela- 
tively long exposure times. We started with a simple ab- 
sorbed power-law fit to each of the spectra of the differ- 
ent sources, separately. The fits were acceptable for 6/10 
sources but residuals at energies less than 2 keV persisted 
in the other 4 sources (NGC 315, NGC 2681, NGC 4278, 
and NGC 5005), suggesting the presence of diffuse hot gas. 

In order to have a better signal to noise ratio and 
photon statistics, we decided to fit the different Chandra 
and/or XMM-Newton spectra of each source simultane- 
ously (the normalizations of the different models in a fit 
were left free between the different EPIC instruments to 
take care of any potential cross calibration uncertainties). 
Simultaneous fitting routines of different observations per- 
formed at different times is done for the first time for any 
sample of LINER sources. For this purpose, whenever a 
source is observed with both Chandra and XMM-Newton, 
we carried out careful imaging analysis of the Chandra ob- 
servations to disentangle the different components (diffuse 
emission, LMXBs, and/or jet emission) that are blended in 



one point-like source in the XMM-Newton extraction region 
of 25"-radius circle (see Appendix B and online Fig. 1-4). 
These components, that we assume non-variable, are in- 
cluded in the simultaneous Chandra/ XMM-Newton fit. We 
do not expec t diffuse hot ga s to vary on a time-scale of 
a few years (Fabbiano 1989), although off- nuclear point- 
like sources and jet X-ray emi ssion could exhibit varia- 
tion on such time-scales (e.g., Harris 2003). We present 
in Appendix B the surrounding medium around the dif- 
ferent LINER Is that are observed with both Chandra and 
XMM-Newton, and we give spectral results to the different 
components. 

We tested multiple models on the data in order to de- 
termine the mechanism responsible for the observed X-ray 
spectra, noticeably: (1) a simple absorbed power-law, (2) 
same as (1) but including a thermal mekal component to 
take account of any diffuse hot gas features in the soft band, 
(3) two power-law components with different photon index 
values with one representing the hard 2-10 keV emission 
and the other representing a possible 0.5-2 keV soft excess 
emissio n commonly seen in the nuclei of luminous galax- 



ies^. g., Porquet et al.|2 004). We investigated spectral vari- 
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Table 4. Best fit parameters to the Chandra snapshot observations of our sample of LINER Is. 



Galaxy Name 


Obs. ID 


N h 


r 


PI Norm, at 1 keV 


kT 


EM< o) 






(10 20 cm" 2 ) 




(10 -5 Photons keV -1 cm -2 s _1 ) 


(keV) 


(10 62 cm" 3 ) 


NGC 266 


1610 


(...) 


1.4 [0.9-1.9] 


2 [1-3] 


(...) 


(...) 


NGC 315 


855 


20 [10-30] 


1.3 [1.1-1.5] 


13 [10-16] 


(...) 


(...) 


NGC 3226 


1616 


74 [47-105] 


1.7 [1.3-2.1] 


22 [15-34] 


(...) 


(...) 


NGC 3718 


3993 


114 [97-132] 


1.5 [1.4-1.7] 


66 [53-81] 


(...) 


(...) 


NGC 4143 


1617 


(...) 


1.9 [1.6-2.1] 


7.2 [6.6-8.3] 


(...) 


(...) 


NGC 4750 


4020 


<31 


1.8 [1.4-2.3] 


5 [3-8] 


(...) 


(...) 


NGC 4772 


3999 


46 [24-53] 


1.69 [1.29-1.74] 


7 [4-11] 


(...) 


(...) 


NGC 5005 


4021 


<39 


2.1 [1.4-2.8] 


5 [3-11] 


0.8 [0.6-1.1] 


3 [1-4] 



Notes. (a) The emission measure (EM) of the meka.1 model, EM= J n e n_j/dV. 

Table 5. Best fit parameters to the LINER Is in our sample observed with a relatively long Chandra and XMM-Newton 
exposure time. 



Galaxy Name 


Obs. ID 


N h 

(10 20 cm- 2 ) 


r 


PI Norm, at 1 keV 
(10" s Photons keV -1 cm- 2 s" 1 


kT 
(keV) 


EM 1 "' 
(10 ,i2 cm- 3 ) 


xl 


d.o.f. 


NGC 315 


4156 


10 [9-13] 


1.5 [1.4-1.6] 


18 [16-21] 


0.55 [0.47-0.59] 


15 [13-17] 


1.13 


465 




0305290201 


L. 


2.1 [1.9-2.2] 


20 [17-25] 


L. 


L. 


NGC 2681 


2060 


<29 


1.5 [1.2-1.8] 


0.6 [0.5-0.8] 


0.67 [0.63-0.70] 


0.4 [0.3-0.5] 


(1.90 


74 




2061 


L. 


L. 


L. 


L. 


L. 


NGC 2787 


4689 
0200250101 


16 [8-24] 
L. 


2.4 [2.1-2.6] 
L. 


3 [2-4] 

4 [3-5] 






1.12 


118 


NGC 3226 


860 


25 [<60] 


1.7 [1.5-2.0] 


13 [9-19] 












0101040301 


89 [82-96] 


1.8 [1.7-1.9] 


25 [23-27] 






0.98 


467 




0400270101 


42 [39-44] 


2.05 [2.0-2.1] 


27 [26-28] 










NGC 3718 


0200430501 
0200431301 


138 [121-155] 
L. 


1.8 [1.7-1.9] 
L. 


57 [48-66] 
47 [40-55] 






0.88 


122 


NGC 3998 


6781 


3 [2-4] 


2.1 [2.0-2.2] 


282 [267-298] 






1.10 


590 




0090020101 


L. 


1.84 [1.82-1.85] 


323 [318-328] 








NGC 4143 


0150010601 


6 [3-9] 


2.2 [2.1-2.3] 


17 [15-19] 






0.98 


118 


NGC 4203 


10535 


(...) 


2.3 [2.2-2.4] 


83 [78-89] 






0.83 


51 


NGC 4278 


4741 


<6.78 


2.1 [2.0-2.3] 


43 [39-47] 


0.62 [0.58-0.66] 


2.6 [2.3-3.0] 








7077 


L. 


2.3 [2.2-2.4] 


18 [17-20] 


L. 


L. 








7078 


L. 


2.3 [2.2-2.5] 


42 [39-46] 


L 


L 


0.93 


310 




7079 


L. 


2.4 [2.3-2.5] 


38 [35-41] 


L. 


L. 








7080 


L. 


2.0 [1.8-2.2] 


11 [10-13] 


L. 


L. 








7081 


L. 


2.1 [2.0-2.3] 


12.5 [11.4-12.9] 


L. 


L. 








0205010101 


3.8[3.1 4.6] 


2.05 [2.03-2.07] 


81 [79-82] 


(...) 


(...) 


1.01 


487 


NGC 5005 


0110930501 


9 [2-18] 


1.7 [1.5-1.8] 


8 [6-9] 


0.64 [0.61-0.67] 


5.2 [4.7-5.7] 


1.2 


107 



Notes. (L.) represents a linked paramter in the fit. ^ a 'The emission measure (EM) of the mekal model, EM= j n e n_ffdV. 



ability in a source by permitting one spectral parameter of 
a given fit to vary independently between different observa- 
tions. We then used the F-test to evaluate the improvement 
in the fit where a 99% probability for an improvement to 
occur by chance is considered valid. The case for NGC 4278 
is already analyzed in Y10 and best fit spectral parameters 
and fluxes are taken from Y10. 

We find that 6/10 sources are best fit with model (1) 
with no additional need for any thermal or two power-law 
components (NGC 2787, NGC 3226, NGC 3718, NGC 3998, 
NGC 4143, and NGC 4203). The rest of the sources were 
best fit with model (2) showing features at low ener- 
gies below 2 keV, indicating emission from diffuse hot 
gas (NGC 315, NGC 2681, NGC 4278, and NGC 5005). 
Model (3) did not improve the quality of the fit in any 
of the cases, giving worse fits at times. We found that 
the intrinsic hydrogen column density varies significantly 



in NGC 3226 decreasing from (8.9 ± 0.7) x 10 21 cm- 2 
to 4.2t°; 2 x 10 21 cm . This is clearly seen in the lower 
panels of Fig. [9] where the soft part of the spectrum dur- 
ing the short XMM-Newton observation is much more ab- 
sorbed compared to the long observation. Additionally, the 
power-law photon index varies in seven sources (NGC 315, 
NGC 3226, NGC 3718, NGC 3998, NGC 4143, NGC 4278, 
and NGC 5005) with the most drastic change being the one 
observed in NGC 315 where T increased from 1.5 ±0.1 dur- 
ing the Chandra observation to 2.1+jj'2 during the XMM- 
Newton one (upper panels of Fig. Mj. This increase is accom- 
panied by a decrease in the 2-10 keV flux from 9.8 x 10~ 13 
to 4.6 x 10~ 13 erg s . This behavior is typical of X-ray 
emission originating in a RIAF structure, which is the ac- 



et al. 


2007) 


tiona 


to L 2 



-to fceV, is inveresly proportional to the photon 
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Table 6. Absorbed fluxes and corrected luminosities derived from the best fit model to our sample of LINER Is and the 
corresponding L 2 -io kev/L E dd- 



Galaxy Name 


Obs. ID 


0.5-2 keV Flux 


2-10 keV Flux 


Corr. 0.5-2 keV Lum. 


Corr. 2-10 keV Lum. 


PI" 


Log(L2-10 keV I LEdd) 








(Logarithmic scale; erg s 1 


cm- 2 ) 




(10 41 


erg s" 1 ) 


% 




NGC 266 


1610 


-13.40 


1 1 K{\ 1 1 Qfll 




-13.10 


-12.80] 


0.23 [0.14 0.29] 


0.59 [0.37 0.74] 


100 


-5.77 


NGC 315 


855 


-12.74 


1 9 78 19 711 
-IZ. 1 O -IZ. 1 ±J 


1199 


-12.03 


-11.95] 


1.56 [1.42 1.71] 


5.42 [4.94 5.94] 


10(1 


-5.42 




4156 


-12.80 


'1909 19 70I 
-IZ.^i -IZ. / J 


-12 04 


-12.05 


-12.02] 


2.42 [2.36 2.53] 


5.06 [4.83 5.18] 


96 


-5.45 




0305290201 


-12.80 


1 9 89 19 78l 

-iz.oz -iz. / 01 


-12 37 


-12.39 


-12.35] 


2.42 [2.36 2.53] 


2.37 [2.26 2.48] 


94 


-5.78 


NGC 2681 


2060 


-13.56 


1 1 ^8 1 1 ^Rl 
-1.1. OO -l.l.ooj 


-13 47 


-13.49 


-13.45] 


0.011 [0.010 0.012] 


0.012 [0.011 0.013] 


73 


-5.80 




2061 


-13.56 


1 1 ^8 11 q^] 
-1.1. OO -lO.-l.Jj 


-13 47 


-13.49 


-13.45] 


0.011 [0.010 0.012] 


0.012 [0.011 0.013] 


73 


-5.80 


NGC 2787 


4689 


-13.30 


1 1 11 11 9fil 
-lO.OO -lO.ZOJ 


-13 28 


-13.32 


-13.25] 


0.0054 [0.0049 0.0058] 


0.004 [0.002 0.007] 


10(1 


-7.70 




0200250101 


-13.34 


' 1 1 'lij 11111 
10.00 IO.OIJ 


-13 33 


-13.37 


-13.29] 


0.0048 [0.0044 0.0052] 


0.0032 [0.0029 0.0034] 


100 


-7.74 


NGC 3226 


1616 


-12.76 


19 81 1 9 71 1 

IZ . 1 IZ . I 1 J 


-12 07 


-12.12 


-12.02] 


0.32 [0.28 0.36] 


0.59 [0.51 0.66] 


100 


-5.38 




860 


-12.76 


1 9 8(1 19 711 
-1Z.OU -IZ. 1 OJ 


-12 32 


-12.36 


-12.29] 


0.19 [0.17 0.20] 


0.33 [0.30 0.35] 


100 


-5.64 




0101040301 


-12.77 


1 9 77 19 7R1 
-iz. / 1 -iz. / OJ 


-12 07 


-12.07 


-12.06] 


0.37 [0.36 0.38] 


0.61 [0.60 0.62] 


100 


-5.37 




0400270101 


-12.56 


19^7 19 


-12 21 


-12.22 


-12.21] 


0.39 [0.38 0.40] 


0.42 [0.41 0.43] 


100 


-5.53 


NGC 3718 


3993 


-12.40 


' 1 9 dl 19 171 

li.ll 1 Z . O M 


-11 50 


-11.52 
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Notes. (a) The power-law component fraction to the 0.5-10 keV corrected luminosity. 



index T (Gu & Cao 2009). Best fit models and residuals 
of the other LINER Is in our sample are shown in online 
Fig. [TO] 

The photon indicies we derived for all of the sources in 
our sample observed with a relatively long exposure time 
varied between 1.5 ± 0.3 and 2.4^Qg with a mean value of 
2.0. Intrinsic column density covered two orders of mag- 
nitude, with Njj varying between ~ 10 20 cm~ 2 for unab- 
sorbed sources, and up to ~ 10 22 cm~ 2 for the only mildly 
absorbed source NGC 3718. The thermal component had 
a temperature mean value of 0.63 keV, consistent with all 
LINER-typ e sources embedded in diffuse emission (Flohic 
eTTaT][2006l ). Table [5] gives the best fit model parameters 
for our sample of LINER Is with a relatively long exposure 
time, and Table [6] gives the 0.5-2 keV and 2-10 keV ob- 
served fluxes and corrected luminosities, as well as the cor- 
responding "Eddington ratios" , L2-10 keV I^Edd- The hard, 
2-10 keV, luminosity spans three orders of magnitude from 



3.2xl0 38 erg s^ 1 to 5.4xl0 41 erg s _1 which resulted in a 
£2-10 kev/LEdd range from 2.0 x 10~ 8 to 2.3 x 10 -5 which 
is at least one to two orders of magnitude smaller than 
L2-10 kev/LEdd seen in lu minous AGN (e.g.: Porquet et al. 
2004] |Nandra et al.|[2007| . 



We looked for any sign of narrow Fe Ka emission line at 
6.4 keV in the EPIC-pn XMM-Newton observations of our 
sample of LINER Is. None of the observations have clear ev- 
idence for the line. Good signal to noise ratio around 6 keV 




Fig. 11. Photon index, T, as a function of the 
£2-10 kev/LEdd 101 our sample of LINER Is. It is clear 
these two quantities are strongly anticorrelated. The solid 
black line designates the least square best fit to a straight 
line. Snapshots, Chandra, and XMM-Newton observations 
are shown in different colours. 



is only found in the spectra of 4 sources that enable the 
estimate of an upper limit: 112 eV for NGC 3718 (obs. ID: 
0200430501), 38 eV for NGC 3226 (obs.ID: 0400270101), 
33 eV for NGC 3998, and 22 eV for NGC 4278. Even though 
a hint for an Fe Ka emission line around 6.4 keV seems to 
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Fig. 12. Top left panel. The positive correlation between the 2-10 keV luminosity and the Eddington ratio. Top right 
panel. The positive correlation between the 2-10 keV luminosity and the BH mass, bottom panels. No correlation is found 
between the BH mass and neither the Eddington ratio (bottom-left) nor T (bottom-right). See text for more details. 



Table 7. Spearman-rank correlation and the Bevington & Robinson ( 2003 1 method to check for the dependence between 
the fit parameters and the intrinsic parameters of our sample of LINER Is. See text for details. 



Correlations 



Coefficient r s Probability (%) Coefficient r Probability (%) 
Spearman-rank Bevington & Robinson (2003 



T — L 2 -10 ketV I ' L>Edd -0.62 

L2-W keV~ L2-IO kev/LEdd 0.61 

1/2-10 fceV - Mbh 0.43 

L2-IO keV / LEdd—MsH -0.41 

T-Mbh 0.33 



>99.9 
>99.9 
99.0 
97.9 
90 



-0.65 
0.65 
0.70 
-0.11 
0.26 



>99.9 
>99.9 
>99.9 

72 

90 



be apparent in the EPIC-pn spectrum of NGC 315, the 
addition of a gaussian line to the best fit model does not 
improve the quality of the fit with an F-test probability of 
60% for an improvement to occur by chance. 



4.3. X-ray correlations 

We looked for any correlations between the X-ray proper- 
ties, mainly the photon index T and the 2-10 keV lumi- 
nosity, and the LINER Is intrinsic parameters, black hole 
mass and the ratio T2-10 kev/LEdd, which could be directly 
linked to the Eddington ratio Lb i/LEd d consider ing that 
Lboi — const, x T2-10 keV (const, sa 16, Ho 2009). We in- 
vestigated the validity of a correlation by fitting the data 
w ith a least sq uare fit to a straight line using the equations 
of York ( 1966 1 and by minimizing the weighted residuals in 



both parameters x and y. We assessed the goodness of the 
fit foll owing the criteria explained in Bevington & Robinson 
( 2003 ) . If a dependent variable y is correlated to a variable 
x with a line slope b, then the reciprocity in fitting x as a 
function of y should lead to a line with a slope b' . Therefore, 
the linear correlation coefficient r, r = Vbl/, varies between 
0, when there is no correlation, to ±1, when there is com- 
plete correlation. This correlation coefficient r cannot be 
used directly to indicate the degree of correlation, instead a 
probability should be calculated that a random sample of N 
uncorrelated data points would yield to a linear-correlation 
coefficient as large as or larger than r. 

We found a strong anticorrelation between the photon 
index T, and L2-10 keV I^Edd-, for our sample. The fit of 
these two quantities with a least square fit to a straight 
line resulted in the following equation: 
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r = (-0.31±0.06)lo 5 (L 2 _ 10 fcey /W) + (0.11±0.40) (3) 

with a linear-correlation coefficient r rs —0.65 and a 
probability greater than 99.99% that L 2 -w kevl^Edd and 
r yield a linear-correlation coefficient > r. Additionally, 
we performed a spearman-rank correlation between T and 
^2-io kev/LEdd and found that these two quantities are 
correlated with a probability greater than 99.99% and a 
spearman-rank correlation coefficient of -0.62. Figure [TT] 
shows the anticorrelation between T and L2-10 kev/LEddi 
and the least square best fit to a straight line. Snapshots, 
Chandra, and XMM-Newton observations are shown in dif- 
ferent colours. 

Using the same criteria as above, we found a positive 
correlation between the hard X-ray luminosity, L2-10 keV, 
and L2-10 keV I LEdd, with a linear correlation coefficient 
r ~ 0.65 and a probability P > 99% that these two quan- 
tities would yield to a linear-correlation coefficient as large 
as or larger than r. Another strong positive correlation we 
found is the increase of the 2-10 keV luminosity with in- 
creasing BH mass with r ss 0.7 and a probability greater 
than 99.9%. We did not find any strong dependence of the 
spectral slope T or L 2 -io kev/^Edd on the BH mass with 
correlation coefficients of 0.26 and -0.11 respectively and a 
probability < 90% that these quantities are correlated. To 
strengthen our conclusions we performed a spearman-rank 
test on the four different correlations and the results are 
shown in Table [7] Our previous results are in agreement 
with the spearman-rank test, except for a weak anticorre- 
lation emerging between £2-10 keV /LEdd and the BH mass 
not seen using our principal correlation test. Fig. |12| shows 
the dependence of the fit parameters to the LINER Is pa- 
rameters. We discuss these results in S 15.41 and S 15.5 



5. Discussion 

In the following discussion, we are comparing X-ray tim- 
ing and spectral results of our sample of LINER Is to 
other results derived on broad samples of LINERs (includ- 
ing type 1 and type 2 LINERs, and transition nucleus) 
and/or low luminosity AGN. We compare our findings, 
X-ray variability and absence of an Fe Ka emission line, 
to type 1 luminous AGN (Seyfert galaxies and quasars). 
We discuss correlations between the fit parameters, T and 
L-2-io keV, and the intrinsic parameters of our LINER 1 
sample, L2-10 kev/LEdd and Mbh, and compare our re- 
sults to luminous AGN and XRBs in order to find out the 
accretion mechanism in LINER Is. 



5.1. X-ray variability 

One of the most important characteristics of AGN is the 
X-ray flux-variability on different time-scales. An anti- 
correlation between the variability amplitude, characterized 
by the normalized excess variance &%xs' with both the 2- 
10 keV luminosity and the black hole mass has been es- 
tablished for a considerable number of AGN ( Nandra et al. 
1997a[ |Turner et al] [19991 |Papadakis||2004l |Q'Neill et al 
2005). Such variability on time-scales of less than a day 
was never detected for LINER sources in the past, with 
observations taken with low spat i al resolution telescopes , 
e.g.:ASCA, (|Komossa et al.||1999[ [Terashima et al.||2002[). 



Ptak et al. ( 1998 ) showed that LINER and low luminos- 
ity AGN sources do not follow the same anticorrelation of 
luminous AGN showing stronger variability with decreas- 
ing luminosity. The authors attributed the non-detection of 
short time-scale variability in LINERs and low luminosity 
AGN to a bigger X-ray emitting region, e.g. RIAF, com- 
pared to luminous AGN. 

Three sources in our sample of LINER Is (NGC 2787, 
NGC 4143, and NGC 4203) show hint of inter-day vari- 
ability with a K-S test probability between 4% and 2% 
that the X-ray emission originates from a constant source. 
Two sources in our sample exhibit significant short time- 
scale variability, both already r eported in the lite rature 
(NGC 4278: Y10 and NGC 3226: |Binder et al.|[2009[ ). 

NGC 4278 exhibited a short time-scale variability (t ~ 
1.5 h) during the XMM-Newton observation where the X- 
ray flux level was highest (compared to the other 6 Chandra 
observations, where no short time-scale variability was de- 
tected). During this 1.5 h time-scale, the flux of NGC 4278 
increased by 10%. Y10 proposed that this variability could 
be the result of a more compact X-ray emission region, e.g. 
an accretion disk truncated at a lower radius during the 
XMM-Newton observation compared to the Chandra ob- 
servations. 

On the other hand, NGC 3226 ([Binder et al 
which is the only source observed for 



2009) 
TOO ks with'XHM- 



Newton, shows significant flux variability during the entire 
observ ation, increasing almost continuously. |Binder et al 



(2009) reported a crfj XS « 0.014, comparable to 0.02 re- 
ported here, which is a variability amplitude similar to 
the one observed in more luminous AGN, but on shorter 



2005) 



Binder et al 



time-scales ( O'Neill et al 
suming a BH mass of 1.4 x 10° 
tio of 2 x 10~ 5 , predicted, using |McHardy et al 
relation, a variability amplitude oF 



(2009), as- 



and an Eddington ra- 
d20Q6b 



2 - a x lO - " on 
a one day time-scale. The discrepancy between the ob- 
ser ved and the pred i cted value of cr NXS 1S the fact that 
the |McHardy et al. (2006) relation is derived for objects 



in a high/soft state. As we show in paragraph 5.4 LINER 



sources, in contrast to luminous Seyfert galaxies, could be 
in a low/hard state similar to XRBs in their low/hard state. 
The L2-10 keV /LEdd derived for the NGC 3226 longest ob- 
servation, where the variability was observed, did not in- 
crease compared to the other three observations, in fact 
it decreased from 4.27xl0~ 6 to 2.95 xlO -6 . Therefore, the 
significant variability seen during the ~ 100 ks observation 
ca nnot be attributed t o an increase in m (as also denied 
by |Binder et aL||2009| ) . |Markowitz k Uttley l ( |2005| showed 
that the PSD break time-scale of the LLAGNNTJC 4258 is 
greater than 4.5 days at > 90% confidence level. The au- 
thors suggested that LLAGN, like XRBs in their low/hard 
state, might have longer break time-scales compared to lu- 
minous AGN and XRBs in the high/soft state. The X-ray 
variability detected in the case of NGC 3226 could be the 
result of a break time-scale of ~ 1 day (~ 10~ 5 Hz) in the 
NGC 3226 PSD, however, this assumption is purely spec- 
ulative and PSD measurement of NGC 3226, which is not 
possible with the present 100 ks observation due to a low 
number of counts for this kind of study, could help confirm 
or discard this idea. Similar to NGC 3226, variability on 
short time-scales (half a day to se veral days) was o bserved 
in the LINER 1 source NGC 3 998 QPian et al.|2010h and in 
the low luminosity AGN M 81 Qlshisa ki et al.|1996||f 
k Makishima||2001 



Pian et al.|2010p . 



lyomoto 
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Fig. 13. Long term X-ray variability of the LINER Is observed more than once. Red ticks represent the error on the 
2-10 keV corrected flux. 



Although all of the normalized excess variance derived 
for our sample consist of upper limits, except for NGC 3226, 
we show in Fig. [8] & 2 nx s as a function of BH mass. Fig. [§] 
rules out any variability on time-scales shorter than 50 ks 
in our sample o f LINER Is, comm only observed in more 
luminous AGN flO'Neill et al.|[2TJ05| , except for NGC 4278 
that shows a 10% dux increase in a ~ 1.5 h time-scale. 

Variability on time-scales of months to years is com- 
mon in our sample. Seven out of nine sources observed for 
more than once show variability on months (NGC 5005 and 
NGC 4278) to years (NGC 315, NGC 3226, NGC 3718, 
NGC 3998, NGC 4143, and NGC 4278) time-scales 8 with a 
5 times flux increase in the most variable source, NGC 4278, 
and a 1.4 times flux increase in the least variable source, 
NGC 4143. The only two sources that do not show vari- 
ability on long time-scales are the ones with a very low 
£2-10 kev/L Ed d- NGC 2681 and NGC 2787 (Table [6}. Both 
sources were observed twice in less than a 5 months pe- 
riod and therefore, variability on years-timescale cannot be 
tested. For illustration purposes, we report in Figure [i"3"| the 
fluxes derived for the sources observed more than once. All 
but NGC 2681 and NGC 2787 exhibit months and/or years 
time-scale variability. 

With the increasing number of X-ray observations of 
LINER and low luminosity AGN sources, more sources are 



All of the sources reported to show variability on years time- 
scale and not on months time-scale lack inter-year observations 
and therefore, monthly variability cannot be tested. 



revealing short time-scale variability (< 1 day). This vari- 
ability is detected with the help of the current generation 
of X-ray telescopes (Chandra, XMM-Newton, and Swift) 
which have a good spatial resolution to isolate the nucleus 
X-ray emission from contaminating sources. Nevertheless, 
and before any firm conclusions are made about variability 
from low accretion rate sources, a homogeneous well sam- 
pled population of LINERs and/or low luminosity AGN, 
with long X-ray observations (~ 100 ks), should be avail- 
able. 



5.2. X-ray spectral shape 



Most recently, |Zhang et al. {[20091) studied the X-ray nu- 
clear activity of a distance limited sample, D<15 Mpc, of 
187 galaxies. The authors found that the X-ray emission 
from ~60% of the elliptical and early type spiral galax- 
ies is consistent with nuclear accreting supermassivc BH. 
They fit the spectra of each source with an absorbed power- 
law and found a photon index L 1.5 — 2 and an intrin- 
sic column density covering almost 4 orders of magnitude 
from 10 20 cm~ 2 to 10 24 cm -2 . The authors found a lu- 



minosity ranging from 



10 38 erg s -1 



to 



10 erg s 



corresponding to a Lo.3_8 kev/^Edd ratio between 10 to 
10~ 8 . Similar r esults were derived in previous work such 
as the work of Gonzalez-Martin et al. (2006 2009) on a 
sample of 82 nearby LINER sources. The authors found 
that 60% of the sample shows a hard unresolved X-ray 
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point-like source in the 4.5-8.0 keV band consistent with 
an AGN assumption. The data were fit mainly by an ab- 
sorbed power-law, with a mean photon index F = 2.1 ± 0.5 
and an intrinsic column density ranging from 10 20 cm -2 to 
10 24 cm -2 , and a thermal component with a mean temper- 
ature kT = 0.5 ± 0.3 keV. 



The results we derive for our homogeneous and optically 
selected sample of LINER Is showing definite detection of 
broad Ha emission are in agreement with all of the results 
derived on broad samples of LINERs and low luminosity 
AGN in the past. X-ray spectra in the whole 0.5-10 keV 
band of the majority of our sample are well fitted with an 
absorbed power-law (9/13). This could mean that the X- 
ray emission in the soft and in the hard band is originating 
from the same region. In the remaining 4 sources (NGC 315, 
NGC 2681, NGC 4278, and NGC 5005), we included a ther- 
mal component to take into account some low energy resid- 
uals, most likely from diffuse hot gas. We found a photon 
index ranging from 1.3±0.2 for the hardest source to 2.41q' 2 
for the softest one with a mean value of 1.9 ± 0.2 and a 

Komossa 



dispersion a = 0.3, similar to values reported in 



(|1999 ) and [Terashima et al. 
density ob 



(2002). The a 



Dsorption 



served in our sample spans two orders 



et al 
column 

of magnitude from 10 20 cm -2 for unabsorbed sources (e.g. 
NGC 3998, NGC 4143) to 10 22 cm" 2 for the mildly ab- 
sorbed source NGC 3718. This is consistent with the fact 
that strong absorption is not expected in objects showing 
broad optical emission line, e.g. LINER Is, as found in lu- 
minous type 1 AGN. T his result was already confirmed by 
Terashima et al. ( |2002 1 who studied the ASCA X-ray obser- 
vations of 21 LINER and 17 low luminosity Seyfert galaxies. 
The authors found a discrepancy in the intrinsic absorbing 
column density between type 1 LINERs, Njj < 10 22 cm~ 2 
(except for NGC 1052 which is not included in our sample, 
see §[2]), and type 2 LINERs absorbed with a column den- 
sity consistent with Njj a few 10 22 cm -2 . The origin of the 
slight excess absorption in some of the sources in our sample 
of LINER Is could be due to either the host galaxy and/or 
some compton thin material intrinsic to the central engine. 
We find a mean value to the temperature of the thermal 
component of 0.63 ± 0.06 keV typical of diffuse hot gas in 
early type galaxies (kT ~ 0.5 — 2.0 keV, Fabbiano 1989). 
Similar kT values were reported in Flofaic et al.| Q2006| ) for 
a sample of 19 LINER sources observed with Chandra. 
The 2-10 keV luminosities of the LINER Is in our sam- 
ple span a range from 3.2xl0 38 erg s _1 to 5.4xl0 41 erg s _1 
which resulted in L2-10 kev/LEdd range from 2.0 x 10~ 8 
to 2.3 x 10~ 5 , in agreement with results derived on broad 
samples of LINERs and low luminosity AGN. These cor- 
responding "Eddington ratios" are at least an order of 
magnitude smaller than those reported f or luminous AGN 
( Porquet et al.|2004| |Nandra et al.||2007 ). 



Although the photon indices measured for our sample 
of LINER Is are similar to those of type 1 Seyfert galax- 
ies (Nandra et al. 1997b I, this does not necessarily mean 
that the X-ray emission in LINER Is is originating from 
an accretion flow s imilar to that of more luminous galax- 
ies (e.g. NGC3998, [Ptak et al.||2004|). Other X-ray timing 



5.3. The absence of an Fe Ka line 

It is now believed that a narrow emission Fe Ka line 
at 6.4 keV is a common feature in the X-ray spectra of 
Seyfert galaxies. The origin of this neutral narrow emission 
line is generally attributed to fluorescence originating from 
parsec-scale dista nces (torus) to d istances closer than the 



broad line region (Shu et al.||2010[). An X -ray Bald win ef- 
fect was discovered by Iwasawa h Taniguchi ( 1993 ) where 



the EW of the Fe line decreases with increasing luminos- 
ity Page et al. (2004) suggested that this Baldwin effect 
observed in X-rays could be the result of a luminosity- 
dependent covering fraction of the putative torus. The in- 
crease in radiation pressure flattens the torus leading to a 
bigger ope ning angle for the torus and, hence, smaller cover- 
ing factor (Konigl & Kartje 1994). This effect was later con- 
firmed on bigger samples of radi o-quiet and radio-loud AGN 



( Nandra et al.|1997c 



Bianchi et al.||2007[ |Chaudhary et al. 
2010 but also sec Jimenez-Bailon et al. 2005) 



LINER Is 

individually studied in the literature with a high signal 
to noise ratio around 6.4 keV do not show any sign of 
Fe Ka emission line with s tringent upper lim its on the 
EW: 25 eV for NGC 3998 flPtak et al.| [20041, 35 eV for 
NGC 3226 ( |Binder et aT]|2009p , and 22 eV for NGC4278 
(Y10). Therefore, it appears that the X-ray Baldwin effect 
do not hold down to very low luminosity AGN. 

We do not detect any significant Fe Ka emission line 
in our sample with upper limits obtained for the sources 
with the highest signal to noise ratio around 6 keV (38 eV: 
NGC 3226, 112 eV: NGC 3718, 33 eV: NGC 3998, and 
22 eV: NGC 4278). If the broad line region is responsible 
for the emission of an Fe Ka line, such as the case for 



the intermediate Seyfert-LINER source NGC 7213 (Lobban 



et al.|2010 Bianchi et al.||2008| ) , we would expect to detect 
it in at least the sources with high signal to noise ratio 
around 6.4 keV, since our sample consists of LINER Is 
showing a definite detection of broad Ha emission. Instead, 
the disappearance of t he torus st ructu re at low Eddington 
ratios as suggested by Ho] (2008) and Zhang et al. (2009) 
could explain the lack of Fe emission lines in our sample, 
assuming that the torus is responsible for the formation 
of the fluorescence Fe line. Interestingly, we find that the 
highest upper limit on the EW of the Fe Ka line is found for 
NGC 3718 which has the highest hydrogen column density 
(N H rj 10 22 cm" 2 ) in our sample of LINER Is. 



5.4. Accretion mode in LINER Is? 

The accretion mechanism responsible for the bulk of energy 
from r adio to X-r ays in LINER sources is still poorly under- 
stood. Ho (2009) recently demonstrated that through local 



and spectral aspects could help shed light on the accretion 
mechanism in LINER Is. 



mass loss from evolved stars and Bondi accretion of hot 
gas, the accretion rate supply needed for the luminosities 
observed in LINERs and other LLAGN is easily attained. 
The author argued that the gas supply present at the center 
of nearby galaxies should generate more active nuclei and 
the luminosity deficit seen in the nearby universe is the 
result of a low radiative efficiency Indeed, radiatively inef- 
ficient accretion flow (RIAF, see Narayan & McClintock 
2008; Quataert|2001 for reviews) models have been applied 
to a growing number of LINERs and LLAGN to explain 
their energy budget an d their spectral energy distribution 



M 81 and NGC 4579 ( Quataert et al 
(iGammie et al.||1999|, NGC 3998 ||Pta' 



1999|), NGC 4258 
k et al. 11200411, and 
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NGC 1097 ( |Nemmen et aLl|2006] ). Wu et al.| ( [2007 ) fit the 
SED of a small sample of eight FR 1 sources with RIAF 



and / or jet m odels and confirmed the prediction of Yuan 



& Cui ( 2005 ) that below a critical value of the Eddington 
ratio (L2-10 keV /LEdd~ 10 -6 ) the X-ray e mission bec omes 



dominated by the jet rather than the RIAF. Ho ( 2009 ) sug- 
gested an even lower value to distinguish between objects in 
the "low" state, where an outer thin disk persists, and those 
in the "quiescent" state containing a pure RIAF. Since we 
are only considering the X-ray properties of this sample in 
our study, we decided to examine the r— £2-10 kev/LEdd 
relation which was shown t o be a good indicator of accre- 
tion rate in luminous AGN ( Shemmer et al.[2006[ ). 

In luminous AGN, a positive correlation is found be- 
tween the hard X-ray power-law slope and the ratio of the 



2-10 keV luminosity to the Eddington luminosity (Wang 
et al.| 2004; Shemmer et al. 2006; Sobolcwska & Papadakis 



2009)7 |Greene fe Ho| ( |2007b| ) showed that the relation holds 
for inter mediate-mass ( 1 ° — 10 6 M ) BHs in active galax- 



ies and Porquet et al. ( 2004 1 showed that the relation 
extends up to more luminous objects when studying a 
sample of 21 low-redshift quasars. A viable explanation 
is that whenever the disk emission increases, the corona, 
the origin of the hard X-ray emission, cools more effi- 
ciently exhibiting a steepening of the hard X-ray spectrum. 
This relation wa s examined for LL AGN (lo cal Seyferts 



and L INERs) by Gu fc Cao (20091 (see also Constantin 



et al. 2009) who found a significant anticorrelation be 



tween the hard X-ray photon index T and the Eddington 
ratio, L bo i/L E dd=30xL 2 ~w kev/LEdd, for the local Seyfert 
galaxies in their sample. However, no strong correlation was 
found when considering only the LINER sources in their 
sample owing most likely to heterogeneous fitting models 
as they have collected their data points from different stud- 
ies. The authors suggested that this anticorrelation found 
in their sample, which is in contrast to the positive cor- 
relation found for more luminous AGN, could signify that 
LLAGN resemble XRBs in th e low/hard state where a sim- 



et al. 


2007 


Wu & Gu 


2008). |Wu & G~ 



not seen in the LINER sample of |Gu fe Cao| ( |2009| This 
strong antico rrelation support the idea suggested by |Gu| 
& Cao (2009) that LLAGN might be similar to XRBs in 
the low/hard state where th e emission is presum ably gen- 
erated in a RIAF structure. Qiao & Liu (2010) predicted 
such an anticorrelation in the low/hard state for their accre- 
tion flow model consisting of an outer-cool optically-thick 
disk and inner-hot optically-thin RIAF within the frame- 



work o f disk and corona w ith mass evaporation ( Liu et al. 
2002|. Qiao & Liu (2010) found that their model can re- 




-3 -2.6 -2 
'°S< L h „/ L ™> 



Fig. 14. Dependence of T as a function of the Eddington 
ratio. The positive and negative correlatio ns represent the 
ones found for a sample of luminous AGNs (Shemmer et al. 
|2008 [) and for our sample of LINER Is respectively. Note 
here that the intercept of our anticorrelation has changed 
since now we are considering L^oi instead of L2-10 keV- The 
crosspoint between the two lines represents a probable tran- 
sitional point from a standard thin accretion disk to a RIAF 
in AGNs. 



dex and the Eddington ratio observed for the X-ray binary 
XTE J1118+480. 



Wu fc Gu ( 2008 1 noted that the transition between the 
two different accretion modes in their sample of XRBs is 
different for different sources but roughly converge to the 
transitional point where T = 1.5 ± 0.1 and log(Lx(Q-5 — 
25 keV)/ LEdd) = —2.1 ± 0.2. Assuming that the accretion 
mode in the high luminosity AGN sample of |Shemmer et ah] 
( 2008 ) is the standard thin disk and the accretion mode in 
our sample of LINER Is is a RIAF, then the transition 
point between the two accretion modes in AGNs woul d be 



14) 



et al. 



sistent with a RIAF whereas a positive correlation could 
mean the existence of the classical thin accretion disk. A 
plausible explanation for the hardening of the spectrum as 
the accretion rate increases in a RIAF context (as seen 
in the r-L2-io keV /LEdd anticorrelation) is the increase of 
the optical depth of the RIAF which leads to an increase 
in the Compton y-parameter resulting in a harder X-ray 
spectrum. 

Our well defined optically-selected sample of LINER Is 
and our homogeneous data analysis techniques allowed 
us to establish a strong anticorrelation between T and 
£2-10 kev/LEdd f° r our sample o f LINER Is (see § 4.3), ([2008]) 



at T « 1.3t£a and log(L bol /L E dd) « -2.6t££ (Fig. 
in good agreement, within the error bars, with the va 
reported for the XRBs and for the sample of |Constantin| 



IK'S 



( 2009 1 . One should keep in mind that these two val- 
ues are affected by several parameters, mainly the BH mass 
and the bolometric l uminosity, that are ca lculated by differ- 
ent me ans between Shemmer et al.| (2008) and this work. 
Indeed, Shemmer et al.| ( |2008p calculated the BH masses 
and the Eddington ratios, L bo i/ Lsddi of their sample of 
luminous AGNs using the vL v (5100k) and FWHM(H/3), 
whereas our calculation of the BH masses of our sample 
are based on the M — a relation and L^oi/LEdd is ca lcu- 
lated using the bolometric correction L bo i = 16 Lx of Ho 



produce the anticorrelation between the X-ray photon in- 



We do not detect any sharp cut below the critical value 
of £2-10 kev/L>Edd of 10 -6 indicating a change in the X- 
ray emissi on, becoming dom inated by synchrotron emission 
from a jet ( Yuan fc Cui|2005[ ). This could be due to selection 
effects where none of our sources, except for NGC 315, com- 
prises strong jets. Moreover, a transition between a "low" 
and a "quiescent" state whenever the accretion rate drops 
below a critical val ue (m s» 10~ 6 ) as suggested most re- 
cently by Ho (2009) could not be tested since, as stressed 
by the author, this kind of analysis should be conducted on 
large samples. 
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5.5. Other correlations and implications 

We found a positive correlation between the hard X-ray lu- 
minosity, £2-10 kev, and L2-10 kev/LEdd- An increase in 
the hard X-ray luminosity with increasing L2-10 kev/LEdd 
is independent of the accretion model and seen in differ- 
ent types of galaxies from transiti on nuclei to high lumi- 
nosity Seyfert nuclei and quasars ( |Ho||2009 Greene & Ho 
2007a). Another strong positive correlation we found for 



our sample of LINER Is is the increase of the 2-10 keV lu- 
minosity with increasing BH mass. A similar trend is seen 
in a sampl e of 112 early- t ype galaxies within a distance 
of 67 Mpc ( |Pellegrinif2010| . Therefore, both the BH mass 
and the £2-10 kev/J^Edd (equivalent to the Eddington ra- 
tio Lboi/LEdd) are driven factors of the hard 2—10 keV 
luminosity. This could be seen in the right upper panel of 
Fig. [12] in § [O] where NGC 2787 with a relatively high 
BH mass of 1.4 x 10 8 M has a low 2-10 keV luminos- 
ity as a consequence of a low L2-10 kev I^Edd (2 x 10 -8 ). 
We did not find any strong dependence of the spectral 
slope r or the Eddington ratio on the BH mass (although 
a spearman-rank correlation showed that £2-10 kev/LEdd 
is weakly anticorrelated to the BH mass). This is in con- 
trast with the anti-correlation between the X-ray spectral 
slope an d the BH mass seen in a sample of low redshift PG 
quasars ( Porquet et al. 2004 ) . It appears that in our sample 
of LINER Is, which have a very low accretion rate, the BH 
mass - here spanning 2 orders of magnitude in \og(MsH) 
space - is not the main driver of the X-ray photon index. 



6. Summary and conclusion 

In the present work, we studied the X-ray properties of a 
sample of nearby LINER galaxies showing definite detec- 
tion of broad Ha emission ( Ho et al.|1997b[ ). Such a sample 
insures the responsibility of accretion into a SMBH for the 
detection of the broad emission lines, guarantees the non- 
existence of large obscuration, and enables X-ray compari- 
son of this class with both XRBs and type 1 AGN. 

Only two sources in our sample exhibit significant hours 
to days time-scale variability. The NGC 4278 flux increased 
by a factor of 10% on a ~ 1.5 hour period, to remain con- 
stant for the rest of the observation (see Y10). On the other 
hand, NGC 3226 shows variability for the whole observation 
of ~ 100 ks, increasing continuously. Three other sources 
show hint of inter-day variability with a K-S test probabil- 
ity between 2% and 4% that the X-ray emission originates 
from a constant source. Short time-scale variability study 
from a homogeneous well sampled population of LINERs 
and/or low luminosity AGN, with long X-ray observations 
(~ 100 ks) should be conducted before any firm conclu- 
sions are made about variability from low accretion rate 
sources. On the other hand, variability on longer (months 
to years) time-scales is common in our sample where 7 out 
of 9 sources exhibit long time-scale variability. The two 
sources not exhibiting variability are the ones with a very 

low L2-IO kev/LEdd- 

The X-ray spectra of our sample of LINER Is are typ- 
ical of all types of LINER sources, fit with an absorbed 
power-law, or a combination of a thermal component and 
an absorbed power-law. We found a photon index for our 
sample between 1.3 ±0.2 for the hardest source and 2.41q' 2 
for the softest one with a mean value of 1.9 ±0.2 and a dis- 
persion o — 0.3. None of the sources in our sample is heavily 



absorbed with NGC 3718 having the highest intrinsic hy 
drogen column density of ~ 10 22 cm~ 2 . The thermal com- 
ponent had a mean temperature kT« 0.6 keV, ty pical of 
other LIN ER sources embedded in diffuse emission (Flohic 
et al. 12006] ) . Our sample spans three orders of magnitude in 
both luminosity and £2-10 keV I^Edd space, ranging from 
10 38 to 10 41 erg s" 1 and from 10~ 8 to 10 -5 respectively, 
which is at least an order of magnitude smaller than the 
Eddington ratios observed in luminous AGN (Seyferts and 
quasars) . 

We do not detect any significant Fe Ka emission line 
at 6.4 keV in the spectra of our sample of LINER Is. We 
obtained upper limits on the Fe line for the four sources 
with the highest signal to noise ratio around 6.4 keV (38 eV: 
NGC3226, 112 eV: NGC 3718, 33 eV: NGC3998, and 22 eV: 
NGC 4278). The lack of a narrow Fe line could be due to the 
disappearance of a torus structure in LLAGN and LINER 



sources ( |Ho|2008|jZhang et al j2009[ ). This implies that the 
X-ray Baldwin effect or the "lwasawa-Tanigushi" effect of 
decreasing Fe Ka EW with increasing 2-10 keV luminosity 
does not extend down to LINER Is. 

Finally, we looked for correlations between the X-ray 
properties and the AGN properties of our LINER 1 sam- 
ple. We found a strong anticorrelation between the power- 
law photon index and the Eddington ratio suggesting that 
LINER Is differ from more luminous Seyfert and quasar 
galaxies that show a positive correlation between the pho- 
ton index and the Eddington ratio. This anticorrelation, 
established for the first time for LINER Is, suggest that 
LINER Is mode of accretion could be similar to that of 
XRBs in their low-hard state, and a RIAF could be re- 
sponsible for the emitting energy from the nucleus. We 
found that the 2-10 keV luminosity in our sample is posi- 
tively correlated to two parameters, the BH mass and the 
Eddington ratio, L2-10 kev/LEdd, confirming the results 
found for broad samples of LINERs and low luminosity 
AGN. On the other hand, it appears that in our sample 
of LINER Is, which have a very low accretion rate, the BH 
mass is not the main driver of the X-ray photon index, as 
the two quantities do not show any strong correlation. 
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Appendix A: Notes on individual sources 



NGC 266. fTcrashima fe Wilson| | |2003| l studied the 2 ks snapshot ob- 
servation made with Chandra and discussed here. They model the 
X-ray spectrum with an absorbed power-law and derive a hydrogen 
column density and a photon index of < 0.82 X 10 22 cm -2 and 1.4, 
respectively. They find a 2-10 keV flux of 1.6 X 10~ 13 erg s~ 1 cm -2 . 
The values of the photon index and the 2-10 keV flux are with a good 
agreement with the value we report here (within the error bars), but 
no additional absorption is added to our mo del. 

NGC 315. Worrall & Birkinsh aw] l |1994[ | first suggested the pres- 
ence of an active galactic nucleus at the center of NGC 315 using 
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ROSAT data. This assumption was later confirmed by |Matsumoto| 
et al. (|2001|) wh en studying a 37 ks ASCA observation (see also 
Terashima et ah| ( |2002^ ). The authors fit the hard 2-10 keV spectrum 
with a power-law and find a photon index of ~ 2 and a luminosity of 
A resolved X-ray jet emission was first reported by 
when studying 



of ft 10 22 cm 
This softening 
3.3 X 10~ 12 erg j3 

NGC 3998. 
observation an> 




2003 ) when studying the Chandra snapshot observation 
e fit to the jet emission with a power-law gave a pho- 
ton index of 2.5±0.7 and a 3.5 X 10 40 erg s" 1 luminosity. The authors 
fit the unresolved core emission with a moderately absorbed power- 
law with an intrinsic hydrogen co lumn density o f ~ 5 X 10 21 cm -2 
and a photon index of 1.4 ± 0.4. |Worrall et al.| ( |2007^ found simi- 
lar results when studying the longer Chandra observation reporting a 
harder core spectrum than th e jet with photon ind icies of ~ 1.6 and 
~ 2.2 respectively. Moreover, |Croston et a l. (2008), after solar flare 
cleaning, studied the only XMM-Newton observation and fit the 60" 
core spectrum with a combination of a mekal and a power-law. Our 
analysis of the two Chandra observations and the XMM-Newton one 
gives similar results to all of the above studies with a simultaneous 
fit to the different extracted spectra where a combination of a mekal 
(kT « 0.5 keV) and a mildly absorbed (N h ft 10 22 cm" 2 ) power-law 
(r between 1.5 and 2) were used. The increase in the power-law slope 
from 1.5 ± 0.1 during the Chandra observation to 2.1^q'J, during the 
XMM-JVewton one is accompanied by a decrease in the 2-10 keV flux 
from 9.8 X 10~ 13 to 4.6 X 10~ 13 erg s _1 . This behavior (see § |5.4| 
is typical of X-ray emission originating in a RIAF structure which is 
believed to be the accretion mechanism responsible for the bulk of 
energy from radio to X-rays in NGC 315 l |Wu et al.|2007[ |. 

NGC 2681 . One of the tw o observations performed with Chandra 
has been already reported in |Satyapal et al.| ( |2005| > . The authors fit 
the 0.5-8 keV spectrum with a combination of a thermal component 
with kT ft 0.7 keV and a power-law with a photon index T « 1.6. No 
in trinsic absorption was required . The same observation was treated 
in |Gonzalez-Martfn et al.| ( |2009| > and same results were derived after 
fitting the spectrum with a mekal and a power-law. iGonzalcz-Martfnl 

I ! ^ I I t . T , ,. - „ f I I P7— • 



nttmg tnc spectrum witn a mekal and a ] 
|et al.| |2009"t derived a 2-10 keV flux of - 



2 x 10~ la 



erg s ± cm 



We fit the spectrum of the two Chandra observations of NGC 2681 
simultaneously with a combination of a thermal component and an 
absorbed power-law. We found similar results to that derived in the 
previous works (T ft: 1.5 and kT ft 0.6 keV) with a 2-10 keV flux of 

3 X 10 -13 erg s -1 cm" 2 . 

NGC 2787. |Ho et al.| | |200T| |, after studying a Chandra snapshot, 
gave this source a class 111 X-ray m orphology, showing hard X -ray 
nucleus embedded in diffuse emission. |Terashima &: W ilson ( 2003]) de- 
rived a 2-10 keV flux of about 3 X 10 4 erg s — 1 cm" a , alter assuming 
a phot on index of 2 and a Galactic absorption. [Gonzalez- Martin et al.| 
J2009J analyzed both Chandra and XMM-Newton long observations 
and fit the Chandra spectrum with a power-law with a rather soft T 
of 2.3 and fit the XMM-Newton spectrum with a combination of two 
power-laws and a thermal component with a power-law absorption 
of ~ 10 22 cm -2 . We found a good fit for both Chandra and XMM- 
Newton spectra simultaneously with a single absorbed power-law and 
found little absorption of ~ 2 X 10 21 cm -2 and a soft power-law pho- 
ton index, V = 2.4, and a 2-10 keV fl ux of 4 X 10~ 14 erg s _1 cm -2 . 

NGC 3226. |George et aT] pOOl"] ) fit the 0.5-10 keV spectrum ex- 
tracted from the long Chandra observat ion with an absorbed ( Nh ~ 
5 X 10 21 cm -2 ) power-law (r ft 1.9). |Terashima fc Wilson| ( |2003| l 
fit the spectrum of the 2.5 ks snapshot Chandra observation with a 
moderately absorbed power-law with Nh ft: 10 22 cm -2 and T ft: 2.2. 
|Gondoin et al.H2004| fit the data of the 35 ks XMM-Newton observa- 
tion with a partial covering absorber to a bremsstrahlung and found 
that the X-ray emitting region, with a temperature kT ft: 0.9 keV 
is 90% covered by an absorber with Njj 



5 X 



10 21 cm~ 2 . |Binder| 
|et al.| ( |2009| l studied the ~100 ks XMM-Newton observation and fit 
the spectrum with a partially covered power-law with T ft 1.9, a 
covering fraction of 90%, and an intrinsic hydrogen column density 
of 10 21 cm -2 . We fit the spectra of both XMM-Newton and the long 
Chandra observations simultaneously with an absorbed power-law and 
found an intrinsic column density varying between the observations 
from ~3x 10 21 cm -2 to ~ 9 X 10 21 cm and a mean photon index 
r = 1.9. 

NGC 3718. |SatyapaI et al.| ( |2005| l studied the snapshot Chandra 
observation and fit the spectrum with an absorbed (Nh ft 10 22 cm -2 ) 
power-law (r ft: 1.5) in excellent agreement with our fit results to 
the same observation. We studied two XMM-Newton observations 
of NGC 3718, being in the field of view of the observations of the 
heavily absorbed Seyfert 2 galaxy UGC 6527. We fit the spectra 
with an absorbed power-law and found an intrinsic hydrogen col- 
umn density similar to the one derived for the Chandra observation 



2 but a somewhat softer power-law with T ft 1.8. 
is accompanied with a 2-10 keV flux decrease from 



|Ptak et al.| p004\ studied the 10 ks XMM-Newton 
d tit the spectrum with a slightly absorbed (Nh ft 



ightly absorbed (Nh 
Same results were found for ob- 



10 20 cm 2 ) power-law (r ft 1.9) 
servations made with BeppoSAX _j|Pellegriniet al,||2000]) and ASCA 
l |Ptak et al.|1999} . |Gon zalez-M artm et al.|(|26f)9} tit the Chandra spec- 
trum with a combination of two absorbed power-laws and a thermal 
component. We fit the XMM-Newton and the Chandra spectra simul- 
taneously with a mildly absorbed (Nh ft 10 20 cm -2 ) power-law and 
found a varying T from 1.8 to 2.1 respectively, occurring with a flux 
decrease from 1 .1 X IP -11 erg s —1 cm" 2 to 6.5 X 10 — 12 erg s — 1 cm -2 . 

NGC 4143. [Terashima fc Wilson| |2003| ) fit the Chandra snapshot 
observation with an absorbed power-law with Nh < 10 21 cm -2 and 
r ft 1.7. We fit the same snapshot observation with a power-law 
without a requirement of an intrinsic absorption and found a similar 
power-law photon index within the error bars, V ft 1.9. We fit the 
XMM-Newton observation with an absorbed power-law and found a 
N H = 6 X 10 20 cm" 2 and a T ft 2.2. 

NG C 4203. A power-law fit to the ASCA spectrum resulted in a 
T ft 1.8 jlyomoto et al.|1998) . [Ho et al.j|200T| | gave NGC 4203 a class I 
X-ray morphology showing dominant X-ray nucleus. We find that the 
40 ks Chandra spectrum is well fitted with a simple power-law affected 
by Galactic absorption with T ft 2.3, softer than the result reported 
for ASCA, most likely due to contamination from X-ray sources in 
the ASCA extraction region of 1'. 

NGC 4278. See Y10. 

NGC 4 750. |Dudik et al.| ( [2005^ and according to the only Chandra 
snapshot observation, gave this source a morphological X-ray type 
II, exhibiting multiple, hard off-nuclear point sources of comparable 
brightness to the nuclear source. We fit the spectrum of this same 
observation with an absorbed power-law with Nh < 3 x 10 21 cm~ 2 
and T = 1.8. 

NGC 4772. We fit the spectrum of the only, Chandra snapshot, 
observation with an absorbed (Nh ~5x 10 21 cm -2 ) power-law (T ft 
1.7). 

NGC 5005. [Terashima et al.| j2002| fit the ASCA spectrum 
with a combination of an absorbed (Nh < 9 x 10 21 c m" 2 ) power- 
law (r ft: 1) and a thermal component (kT ft: 0.8). |Dudik et al.| 
(2005) gave NGC 5005 a morphological X-ray type III, showing a 
hard nuclear point source embedded in diffuse emission. | Gonzalez- 1 
|Martin et al.| l |2009[ ) fit the XMM-Newton spectrum with a combi- 
nation of a thermal component with kT ft: 0.3 keV and an absorbed 
(Nh ~ 6 X 10 21 cm -2 ) power-law (r ft 1.5). We find a hotter thermal 
component when fitting the same data set with kT ft: 0.6 keV and a 
mildly absorbed power-law with Nh ft 10 21 cm -2 and T ft: 1.7. 



Appendix B: Surrounding sources of the centers of 
galaxies observed with Chandra 

In this Appendix, we report the spectral analysis of resolved and/or 
unresolved off-nucleus sources detected with the Chandra telescope, 
but blended w ithin the central LINER in the XMM-Newton extrac- 
tion region (§ |3.2| . Fig. 1-4 show the surrounding sources of these 
LINER Is detected with Chandra within a 25"-radius circle. The sur- 
rounding medium of NGC 4278 is already reported in Y10. 

NGC 315. This source is the only source in our sample that shows 
a resolved X-ray jet. We extracted from the longer Chandra observa- 
tion the spectrum of the jet from an ellipse with a semi-major axis 
of about 11.3" and semi-minor axis of 5.6". The base of the ellipse 
extends down to the 1.1 X 99% PSF of the central source. We fit the 
spectrum with a combination of an absorbed power-law and a thermal 
mekal component and found a good fit with a reduced \ 2 °f 0-9 for 42 
d.o.f. We find a hydrogen column density upper limit of 2 X 10 21 cm~ 2 



and a photon index Y = 2.0_ 



The thermal component had a tem- 



perature of 0.6^Q g keV. We find a corrected 0.5-10 keV flux for the 
jet emission of (1.2 ± 0.1) X 10 — 13 erg s — 1 cm -2 , which corresponds 
to a 0.5-10 keV luminosity of 6 X 10 40 erg s — 1 , corresponding to 8% 
of the nuclear flux. The power-law emission contributes to 95% to the 
total emission of the jet. For the diffuse emission, we extracted the 
spectrum from an annulus with inner circle delimited by 1.1 X 99% 
PSF of the central source and an outer radius of 25" excluding the 
jet extraction region. The same model fit to the jet gave a good fit 
with a reduced x 2 of 1.3 for 61 d.o.f. We find an intrinsic absorption 



to the powcrlaw component Nh 



cm and a photon 
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index T = 1.7 ± 0.8, possibly representing emission from unresolved 
X-ray binaries. The thermal component has a temperature similar to 
the one derived for the jet emission with kT = 0.6 ± 0.2 keV. We 
found a corrected 0.5-10 keV flux of 2.9±o' 2 x 10" 13 erg s" 1 cm -2 
with the power-law contributing only to 30% of the total emission. 
This corresponds to a luminosity of ~ 10 41 erg s — 1 which is ~ 14% 
of the total core luminosity. 

NGC 2787. An X-ray source south-east of the nucleus of 
NGC 2787 at a distance less than 10" is present. We fit the spectrum 
of this source with an absorbed power-law. We used the same red- 
shift and Galactic absorption as for the NGC 2787 nucleus, assuming 
that this X-ray source is located in NGC 2787 and not a background 
quasar. The fit is acceptable with a reduced x 2 °f 1-0 for 19 d.o.f. We 
found a power-law photon index T = 1.5^g';j, typical of X-ray binaries 
in nearby galaxies jlrwin et al.| 2003 Fabbiano 2006), and an upper 
limit on the intrinsic hydrogen column density ot 2 X 10 21 cm -2 . We 
derived a 0.5-10 keV corrected flux of 7 ± 1 X 10 — 14 erg s — 1 cm -2 , 
which resulted in a luminosity of, adapting the NGC 2787 distance 
of 7.48 Mpc, 5 ± 1 X 10 38 erg s — 1 . This luminosity is close to the 
NGC 2787 core luminosity of ~ 9 X 10 38 erg s — 1 . Such a source could 
be a luminous low mas s X-ray binary ( LMXB) similar to some seen 
in early type galaxies ifFibbiano 20061. The rest of the medium in 
an annulus of outer radius 25" around NGC 2787 is formed by six 
other X-ray sources, much fainter than the closest one to the cen- 
ter. We could not perform spectral analysis on the different sources 
aside, so we fit the spectrum of the six sources simultaneously with a 
power-law only affected by Galactic absorption. We found a photon 
index of ~2 for the six sources and a total corrected 0.5-10 keV flux 
of 8±2 x 10~ 15 erg s 1 cm 2 , resulting in a 0.5-10 keV luminosity of 

5^2 x 10 37 erg s —1 , corresponding to 5% of the nucleus luminosity, 
when adapting the NGC 2787 distance of 7.48 Mpc. 

NGC 3226. Two X-ray sources are within a ~12" distance from 
the nucleus of NGC 3226 (source 1: CXOU J102334.1+195347 , source 
2: CXOU J1 02326.7+195407). Both sources were reported in [George] 
|et al.|fl200i) . Based on the hardness ratio between the counts in the 
0.3-2 keV band and the counts in the 2-10 keV band, the authors 
estimated the sources to have a flux of a few 10 -14 erg s — 1 cm -2 , and 
therefore a luminosity between a few times 10 38 erg s _1 to a few times 
10 39 erg s — 1 . We fit the spectrum of both sources with an absorbed 
power-law, using the cash statistic due to a low number of counts. We 
find that source 1 and source 2 have an intrinsic absorption upper limit 
of 6 X 10 22 cm -2 and 5 X 10 21 cm -2 , respectively. The photon index 
we found for source 2 is typical, within the error bars, of accreting 
objects with a T = 1.2^2'^. On the other hand, we found a harder 
spectrum for source 1 with a T = 0.2 ± 1.3. The corrected 0.5-10 keV 
flux we derive for source 1 and source 2 are ~ 4 X 10 -14 erg s" 1 cm" 2 
and ~ 2 X 10 — 14 erg s — 1 cm -2 , respectively. This implies, assuming 
the distance of NGC 3226 to both sources, a corrected 0.5-10 keV 
luminosity of 3 X 10 39 erg s — 1 and 2 X 10 39 erg s _1 for source 1 and 
source 2 respectively. Both luminosities are well beyond the luminosity 
of a typical neutron star LMXB of ~ 3 X 10 38 erg s _1 , and hence 
could be BHs greater than or equal to a few solar masses, very young 
supernovae, or microquasars (Geor ge et al.|2001^ . 

NGC 3998. Only one source is detected within a 25" circle around 
NGC 3998 in the Chandra image. We fit the source with an ab- 
sorbed power-law, using the cash statistic. We found an upper limit 
on the intrinsic column density of 10 21 cm -1 and a power-law pho- 
ton index of 1.4^g' 9 . We derived a corrected 0.5-10 keV flux of 
3 i 1 X 10 — 14 erg s — 1 cm -2 , which corresponds to a luminosity of 
7±lxl0 38 ergs" 1 , adapting the NGC 3998 distance of 14.1 Mpc. This 
corresponds to only 0.2% of the total core luminosity of NGC 3998 
and match the luminosity of XRBs in nearby galaxies. 
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Fig. 1. Chandra image of the central 25" of NGC 315. Jet Fig. 4. Chandra image of the central 25" of NGC 3998. 
spectrum is extracted from an ellipse with a semi-major Only one source, much fainter than the core, is present in 
axis of 11.3" and a semi-minor axis of 5.6 ". Core emission the field. See Appendix B for more details. The horizontal 
is extracted from a circle centered on the source with a bright line corresponds to the readout streak events, 
radius comprising 99% of the PSF (—2.7"). The rest of the 
medium inside a 25" circle is considered diffuse emission. 
See Appendix B for more details. 




Fig. 2. Chandra image of the central 25" of NGC 2787. 
A point-like source south-east of the central LINER is 
present with a luminosity comparable to the core luminos- 
ity. Another six point-like sources, marked in white, arc 
present in the field. See Appendix B for more details. 




Fig. 3. Chandra image of the central 25" of NGC 3226. 
Two sources are present in the field marked source 1 and 
source 2. See Appendix B for more details. 
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Fig. 6. Light curves and hardness ratios of the LINER Is observed with XMM-Newton, all binned with a 1 ks time 
bin-size. 
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Fig. 7. Light curves and hardness ratios of the LINER Is observed with Chandra with a long exposure time, all binned 
with a 1 ks time bin-size. 
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Fig. 10. Data and best fit model of the spectra of the LINER Is in our sample with a relatively long exposure time. 
Residuals of every fit are given in terms of sigma. 



